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Abstract 

We investigate magnetic properties and statistical effects in ID strongly repulsive two-component 

fermions and in a ID mixture of strongly repulsive polarized fermions and bosons. Universality 

in the characteristics of phase transitions, magnetization and susceptibility in the presence of an 

external magnetic field H are analyzed from the exact thermodynamic Bethe ansatz solution. We 

show explicitly that polarized fermions with a repulsive interaction have antiferromagnetic behavior 

at zero temperature. A universality class of linear field-dependent magnetization persists for weak 

and finite strong interaction. The system is fully polarized when the external field exceeds the 

critical value ~ ^Ep, where Ep is the Fermi energy and 7 is the dimensionless interaction 

strength. In contrast, the mixture of polarized fermions and bosons in an external field exhibits 

square-root field-dependent magnetization in the vicinities of H = and the critical value H = 
H M ra \& Ep We find that 

a pure boson phase occurs in the absence of the external field, fully- 
polarized fermions and bosons coexist for < H < , and a fully-polarized fermion phase 
occurs for H > . This phase diagram for the Bose-Fermi mixture is reminiscent of weakly 
attractive fermions with population imbalance, where the interacting fermions with opposite spins 
form singlet pairs. 

PACS numbers: 03.75.Ss, 05.30.Fk, 71.10.Pm 
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I. INTRODUCTION 



The recent success in experimentally realizing degenerate quantum gases in low di- 
mensions 

a a a a a 

has revived interest in one- dimensional (ID) integrable mod- 
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131 ] . The ID atomic 



els of interacting fermions and bosons aaaa 

gases with internal degrees of freedom are tunable interacting many-body systems featur- 
ing novel strong correlation effects and subtle quantum phase transitions. Exotic quan- 
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rases in ID two-component attractive 
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ermions have received considerable interest 
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28| due to the experimental observation 
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3jJ. For certain regimes, the two- 



of fermionic superfluidity and phase transitions 
component Fermi gas with population imbalance can be viewed as a strongly interacting 



Bose-Fermi mixture 
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331 ] . However, subtle differences between the bosonic dimer and 



pure bosons have been observed experimentally 3J, [35|. For repulsive interaction, the two 



14 



15, 



36 



4fJ, 
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39] . In contrast 



421 ] possesses novel 



component Fermi gas exhibits antiferromagnetic behavior 
to the two-component Fermi gas, the ID spinor Bose gases |10. 
ferromagnetic properties and collective dynamics of spin waves at temperatures below the 
degenerate temperature. The subtlety in ID quantum many-body physics [43J is a result of 
the dynamical interaction between identical particles and their statistics. 

On the other hand, the recent success in tuning interspecies Feshbach resonances in 
Bose-Fermi mixtures of ultracold atoms 
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44|, such as 6 Li+ 7 Li, 6 Li+ 23 Na, 40 K+ 87 Rb, 



6 Li+ 87 Rb, opens up a new gateway for exploring many-body physics, including superflu- 
ids and Mott insulators, spin and charge density waves, phase separation, the BCS-BEC 
crossover etc. On the theoretical side, Bose-Fermi mixtures have been studied through var 
ious techniques like the mean-field approximation 



Quantum 



ansatz 



51 



^ Luttinger liquid formalism [46j, |47| . 
Q, 5(3] > exac t solutions using the Bethe 
521 ] and other methods 53] . However, the role of quantum fluctuations is en- 



Vlonte Carlo [48j , bosonization techniques 



hanced in ID compared to the three-dimensional case to the extent that traditional mean- 
field theories fail for strong interaction in ID. The exact Bethe ansatz solution of ID many- 
body systems provides more reliable physics than the mean-field theory. 

In this paper, we investigate external field-dependent magnetic properties and statisti- 
cal effects of ID two-component fermions with repulsion and a ID mixture of polarized 
fermions and bosons by means of their exact thermodynamic Bethe ansatz solution. The 
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antiferromagnetic groundstate properties and critical behavior of two-component repulsive 
fermions are studied in detail. The universality class of linear field-dependent behavior of 
magnetization is predicted for weak and finitely strong coupling regimes. However, in the 
Tonks-Girardeau limit, a van Hove singularity in the susceptibility occurs as the external 
field approaches a critical field. In this limit, the system becomes a paramagnet, i.e., a 
system with free spins. However, the existence of bosons in the Bose-Fermi-Fermi mixture 
destroys the antiferromagnetic ordering in the two-component interacting fermions. In the 
absence of the external field, the state of pure bosons is among the groundstates. When the 
external field is turned on, a second order phase transition from a pure boson phase into a 
mixed boson-fermion phase occurs. A fully-polarized Fermi liquid occurs when the external 
field exceeds a critical value. We calculate the explicit details of these transitions. 

This paper is set out as follows. In section [Til we recall the Bethe ansatz solution of the 
ID integrable model of mixed bosons and fermions. In section IIII} we study the antifer- 
romagnetic behavior of the ID strongly repulsive interacting Fermi gas with polarization. 
The phase diagram of the model is presented. The quantum phase transitions and magnetic 
properties of the Bose-Fermi mixture are studied in section [IVl Section [V] is devoted to 
concluding remarks and a brief discussion. 



II. THE HAMILTONIAN 



We consider a 5-function interacting system of TV bosons and fermions with equal mass 
m and internal degrees of freedom, or a mixture of two-component fermions and spinless 
bosons, constrained by periodic boundary conditions to a line of length L in an external 
magnetic field H. The Hamiltonian 0, 55] is 

i=l % l<i<j<N 

-Ih(m 1 -m 2 ). (i) 

The quantum numbers Mi and M 2 are the numbers of fermions with spin-up and spin-down, 
respectively. Under exchange of spatial and internal spin coordinates between two particles 
the wavefunctions of the Hamiltonian ([1]) are symmetric for bosons or for fermions with 
opposite hyperfine states and antisymmetric for fermions with the same hyperfine states. 
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For maintaining the integrability of the model, we consider the same mass for bosons 
and fermions and the same interaction strength gm between bosons, between bosons and 
fermions and between fermions with opposite hyperfine states. The interaction is attractive 
for gm < and repulsive for gm > 0. There is no 5-type interaction between fermions with 
the same hyperfine states. Although these conditions appear somewhat restrictive, the model 



sti 



22, 



1 captures the essential physics 
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5l| relevant to the current experimental observations 
351 ] . The interspecies interaction can be tuned from strongly attractive (gm — > — oo) 



to strongly repulsive (gm — > +oo) via Feshbach resonances and optical confinement. The 
coupling constant gm can be written in terms of the scattering strength c = 2 /am as 
fi'iD — fi 2 c/m. In principle an effective ID scattering length am can be expressed through 
the 3D scattering length for the bosons and fermions confined in a ID geometry. We use 
the dimensionless coupling constant 7 = mgm/(fo 2 n) for physical analysis. Here n = N/L 
is the linear density. 

This model was solved in 1963 by Lieb and Liniger for the special case of spinless bosons 
56]. Experimental realization of ID interacting Bose gases [l|, y, [3], 4] has stimulated further 
interest in various integrable models, including interacting two-component fermions solved 



by Yang and Gaudin 



57|, 



581 ]. the Bose- Fermi mixture solved by Lai and Yang 55], and 



multi-component bosons and fermions solved by Sutherland |54|. 

The Bethe ansatz equations for the Bose-Fermi mixture ([1]) with an irreducible represen- 
tation [2 + M b , 2 Ah ~\ l M i" M 2] are [55] 



M 



N \ 

n 



exp(ifcjL) 
— kj + ic' 



j-j kj — X a + id 



kn 

Q = l 3 

M 



n 

13=1 



A Q — A/3 + ic 
A« — A/? — ic 



x 



n 



ic 



b=i 



\ a - A b + id 



(2) 



M 

n 

k=l 



A b - A fc - ic ; 
A b - X k + id 



1. 



We define M = M 2 + M b , where M b is the number of bosons. In these equations kj, 
with j = 1, . . . , N, are the quasimomenta of the particles and X a , with a = 1, . . . , M, are 
parameters for the fermions with spin-down and bosons. A b with b = 1, . . . , M b are the 
parameters for the bosons. In the Bethe ansatz process [55] . the conjugate representation 



[Mi, M2, 1 2+Mfc ] was used to solve matrix eigenvalue equations associated with the irreducible 



representation [2 + Mb, 2 M2_1 , i A/ i~ A/ 2] following Sutherland's paper 54j. 

The Bethe ansatz solution of the Bose- Fermi mixture was considered with renewec 
terest 
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m- 



52| due to the experimental observation of new quantum phases 34], |35 |. 
The groundstate properties, correlation functions and harmonic trapping for the mixture of 
fully-polarized fermions and bosons were studied in Ref. [5l|. The model of arbitrarily po- 



larized fermions and spinless bosons was investigated in the weak coupling regime 



52j. The 



Bethe ansatz equations (j2J) also contain those for two-component fermions [57| when M& = 0. 
However, the magnetic properties and phase transitions for two-component fermions in the 
repulsive regime and the mixture of polarized fermions and bosons have not been compre- 
hensively investigated. We turn now to the investigation of their magnetic properties and 
phase transitons in the repulsive regime. 

III. TWO-COMPONENT FERMIONS 

It is well known that the groundstate for ID interacting fermions is antiferromagnetic. 
This is proved to be the case by the Lieb-Mattis theorem [5^. They showed that for a 
system with total spin S' > S, the lowest energy belonging to that system is also greater, 
E(S') > E(S). Hence when there is no external field, the system is not polarized. Another 
way of showing this is from the Pauli exclusion principle. When there is no interaction at 
all, the groundstate can only have one pair of fermions (spin-up and spin-down) having a 
particular pseudomomentum k. When a repulsive interaction 7 is "switched" on, the pair 
will split so that no two fermions can have the same value for k. However, the sequence 
of spin-up and spin-down fermions is still unaltered, which shows that they remain anti- 
ferromagnetic even when 7 > 0. The magnetization of the two-component fermion system 
in the weak coupling limit, 7 < 1, can be directly obtained from the groundstate energy 



derived from the discrete BA equations 



37] 



£ ~ £ n 2(l _ p2) + ^ + ^ln 3 P 2 (3) 
L 2 1 12 4 y 1 

where in an obvious notation the polarization P = (N^ — N^/N. Defining the magnetization 
per particle as m z = nP/2, the magnetization has a linear field- dependent form given by 

H « 2n (tt 2 - 2 7 ) m z . (4) 

5 



When the external field exceeds the critical value ~ n 2 (ir 2 — 27) a fully-polarized phase 
occurs. The magnetic susceptibility is % ~ l/2n(7r 2 — 27) for finite H . However, from field 
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theory, \ ~ V(2 7r ^s) at if = 0. Here the spin velocity is v s ~ nir(l — 7/7r 2 ) 

We turn now to the magnetic properties of the ID exactly solved model of two-component 
fermions with arbitrary polarization in the strong coupling regime 7 ^> 1. In the thermo- 
dynamic limit, L, N — > 00 with N/L finite, the Bethe ansatz equations can be written in 



terms of the spin and charge densities in the form 



57 



p(k) + p h (k) 



1 1 

2^ + 2^ 



B 



ca{\) 



B c 2 /4 + (k - A) 2 
cp(k) 



dX 



(5) 



c 2 /4 + (A - k) 2 
B 2ca(X) 



dk 
dX'. 



(6) 



2vr J_ b c 2 + (A - A') 2 

The integration limits Q and B are determined by the total number of particles N/L = 
J-q p(k)dk and the number of spin-down fermions M 2 /L = J B B cr(\)d\. For the groundstate 
p h (k) = (T h (X) = 0, i.e. there are no holes for both charge and spin rapidities. 

The groundstate properties can be determined from the TBA equations in the limit 
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621 ] . In terms of the dressed energies e(k) := Tln[p h (k)/ p(k)] for the charge and 



T -> 

4>(X) := T\n[a h (X)/a(X)} for the spin degrees of freedom, these equations are 



e(Jfe) = k 

0(A) = fc(A) 
where the driving term is given by 



B 



c0-(A) 



M 2 2tt J_ b c 2 /4 + (k - A) 2 



rfA 



1 

2^ 



2c0-(A') 



^i(A) 



— 



B c 2 + (X- A') 2 
ce-(ife) 



(7) 
(8) 

(9) 



2tt 7_ q c 2 /4 + (A - k) 2 

The negative part of the dressed energy, e(k) < for \k\ < or 0(A) < for |A| < 
B, corresponds to occupied states, with the positive part of e(k) and 0(A) corresponding 
to unoccupied states. We clearly see from the TBA equations (j7J) and (jHJ) that the spin 
interaction is antiferromagnetic with an effective spin-spin exchange interaction depending 
on the energy of the system. For H = 0, the driving term gives rise to an asymptotic 
condition 0(oo) = below the A-axis. This leads to a maximum 0~(A) which gives the 
lowest energy state of the system. For strong coupling 7 > 1, the driving terms in the 



second TBA equation © become (71(A) ~ H — cPq/(c 2 /A + A 2 ) where we ignore 0(l/c 2 ) 
contributions. Here Pq = J^e~ (k)dk is the pressure. The driving term gi clearly 
indicates that the form of equation (jBJ) indicates an effective antiferromagnetic spin-spin 
exchange interaction with an effective coupling constant J = —2Pq/c < in contrast to 
the effective ferromagnetic coupling in spinor Bose gas Q|, for which J = 2Pq/c. For the 
groundstate we find that (/>o(\) = — c cos ^(° A i c \ ■ On solving the TBA equation (jSj), the leading 
behavior for the energy per unit length and the chemical potential is given by 

1 % 9 ( 41n2\ 
F(0)/L = -n 3 vr 2 1 (10) 



3 V 7 
6 In 

^7 



^nV(l-^l (11) 



which coincide with previous results 



A. Finite external field H 



l4. 



When the external field is applied, the spin-up states are energetically favoured due to 
Zeeman splitting. As the magnetic field increases, the Fermi sea of the spin sector (f>(X) is 
raised. For small field, H <C 1, the second TBA equation (JSJ) becomes 

0(A) = \H 2rf 



with the function 



2 2ccosh(7rA/c) 

/oo 
G(X - X')(f) + (X')dX' (12) 
■00 

G(X) = — / r^e-^du. (13) 

y J 2tt 1 + el— I V ; 



Making a change of variables A — > A + B and defining y(X) = 0(A + B), equation (|T2 
becomes 

y(A) = oH~- 7—^ + / G(A - k)y(k)dk 

2 2c cosh(7r(A + f>)/c) Jo 



/"OO 

+ / G(A + 25 + k)y{k)dk. (14) 

JO 
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Observing that G(X+2B+k) — ► e 2nB / c for B > 1, we can take the expansion y = yi+y 2 +- 
with respect to the order of e - 2nmB / c where m 

1 



1,2, 



It follows that 



yi(A) + y 2 (X) + ... = -H- 2TtPc- 1 e-^/ c e-* B ' c 



+ 
+ 

Collecting terms of the same order then gives 

1 



roc 

/ G(X + 2B + k)(y 1 (k)+y 2 (k) + ...)dk 
Jo 

G(X-k)( yi {k)+y 2 {k) + ...)dk. 



Vi(A) 

y 2 (A) 



# - 27rP(rV" A/c e-" B/c + / G(X - k) Vl {k)dk 

/»oo /*oo 

/ G(A + 25 + k) yi (k)dk+ / G{X - k)y 2 {k)dk 
Jo Jo 



(15) 

(16) 
(17) 



and so on. To a sufficient approximation, we only need to consider the equation for yi(X). In 
considering the leading contribution in powers of the exponential G(X + 2B + k) — > q- 27tB / c 
for B 3> 1, equation (fl2|) reduces to the standard Wiener- Hopf equation (for the Wiener- 



Hopf technique, see, e.g., [62 
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energy, pressure and chemical potential as 



F 
P 



-n 2 n 3 I I 



67|). After tedious calculation we obtain the free 
3H 2 ~f A 12 In 2 s 



41n2 

7 



1 + 



7r 2 n 3 



61n2\ 9# 2 7 , 

1 I + "a 1 + 



7 



^vr 4 n 



., , , 161n2\ 3F 2 7 , , 
^ sa 7r 2 n 2 [ 1 ) + I 1 + 



37 



4tt% 2 



7 

104 In 2 

97 
34 In 2" 

37 



(18) 
(19) 
(20) 



Moreover, the magnetization m z (H) and susceptibility x{H) follow in the small field regime 
using the formula m z = —dF(H)/dH and x{H) = dm z /dH. Hence 



m z {H) 
X(H) 



3# 7 




12 In 2 


47r 4 n 2 


7 


37 




12 In 2 


4ttV 


7 



(21) 
(22) 



In the small field regime, the higher-order logarithmic correction to the susceptibility can 
be calculated in the same fashion by solving to next order in the Wiener-Hopf equations. 
Thus finite logarithmic singularities are expected in the susceptibility 62j, I64J. As n — > the 
susceptibility diverges as a consequence of the van Hove singularity of the empty band. The 
linear field-dependent behavior of the magnetization is maintained for finitely strong and 
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weak repulsion, see figured! However, for the Tonks- Girardeau limit, 7 — > 00, caution should 
be paid to the critical behavior of the polarized fermions. The spin velocity is proportional 
to the inverse of the susceptibility. Therefore, in this limit the spin velocity v s tends to 
zero such that the spin propagation is almost frozen. As a consequence the statistical effect 
is fully suppressed and the Tonks-Girardeau Fermi gas becomes a paramagnet or say free 
spins. Therefore, we expect an infinitely divergent susceptibility to occur for the system 
in the limit 7 — > 00. In the next part, we shall search for evidence of the paramagnetic 
behavior in the Tonks-Girardeau Fermi gas with population imbalance. 

If the external magnetic field is sufficiently strong, the system can be fully-polarized. 
From the the TBA equations and (151), we find that if the field H > = 4P /c, then 
4>(0) > 0. The critical value indicates a phase transition from a partially-polarized 
ferromagnetic phase to a fully-polarized ferromagnetic phase. For the limit H — > H^, the 
system is almost fully polarized, corresponding to a very small Fermi sea for the spin sector 
of the TBA equation. The critical behavior may be directly derived from the TBA equations 
(jZJ) and (jSJ) or the Bethe ansatz equations (jSJ). In the vicinity of H^, we can approximate 
the integral for the spin part by taking the area under the curve as a rectangle as usually 



done in spin chains and ladders 



651 ] . In this way 



p c 2Bc(p-(0) 

with 0(0) ~ 1+23/^0 • Further, (f)(B) = provides the value of the integration boundary, 
with 



On the other hand, from the Bethe ansatz equations (JBJ), we find 



1 



B = -C7r(l - Vl - 2a) (25) 

if a = M/N 1. Here M = M2 is the number of spin-down fermions. Thus we have an 
explicit relation 

a = -y/(H?-H)/(5H?-H) (26) 

7T 

between the numer of spin-down fermions and the external field, where = 4P/c. 

Substituting the expression (|23|) into the first equation of the TBA ((SJ) and also using the 
conditions e(Q) = and dP/dp, = n, we obtain the leading terms for the pressure and the 



free energy 



The groundstate energy 

L 3 " " V' 7 

can be obtained from the relation E/L = F(H)/L + Hn(l — 2a)/2, which is consistent 



Ion f 8a\ 
« ^2 i _ (29) 



with the BA result 



37|] . Here recall that a = M/N. The critical field iff « 4P/c « ^E F 



where Ep = |n 2 7r 2 is the Fermi energy. These explicit results indicate that for 7 — > 00, the 
system becomes a paramagnet and statistical interaction is suppressed. It also follows that 
the magnetization m z (H) and susceptibility x{H) in the limit H — > are given by 

m ^~2 n(5H-H e y/* n(5H-H c )V 2 [M) 



X(H) 



3n 



ir(5H - H c ) l / 2 {H - H c ) 1 / 2 
lhn(H c - Hf/ 2 



(31) 



n (5H - H c ) 3 / 2 ' 

The expression for the susceptibility is seen to be divergent when H — > . A quantum 
phase transition occurs as the external field is greater than the critical field H^, see the 
phase diagram in figure El The square-root behavior of the magnetization is a consequence 
of the van Hove singularity in one dimension. 

The analytic results for the magnetic properties show that spin-spin interaction for system 
of the polarized fermions with strong coupling 7 > 1 can be effectively described by the 
isotropic spin- 1/2 Heisenberg chain with a weak antiferromagnetic coupling J « —AEp/^f. 
For the isotropic spin-1/2 Heisenberg chain 62], the susceptibility at H = is given by 
X ~ 1/Jtt 2 which coincides with the result (1221) . Furthermore, we can find the spin velocity 
at H = 0, 

2A <l - ^) . (32) 



37 V 7 

The antiferromagnetic spin-spin exchange interaction embeded in the ID interaction po- 
larized fermions may provide insight into understanding the spin segregation in a trapped 
Fermi gas 
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IV. BOSE-FERMI MIXTURE 



The earliest study of the groundstate of the ID mixture of bosons and polarized fermions 



dates back to Lai and Yang in 1971 



55| who showed by numerically solving a set of coupled 



integral equations that the energy of the system is a monotonic decreasing function with 
respect to the ratio of bosons and the total number of particles in the system. This implied 
that the groundstate is occupied only by bosons. The same signature was found recently by 



other groups 50, 



51 



explicitly in Ref. 



521 ] . The groundstate energy for the weak coupling regime was given 



52j. In addition, study of the quasimomentum distributions |68| showed 



that the quasimomentum distribution is more compressed when there are no fermions. As 
the number of fermions increases, the quasimomentum distribution becomes more spread 
out due to the Fermi pressure. 

It is thus natural to ask what the true physical groundstate is in the presence of an 
external magnetic field. What kind of magnetism does the model exhibit? As far as we 
understand, the magnetic properties of the mixture of bosons and polarzied fermions have 
not been studied yet. It is of interest to know how the addition of bosons influences the 
antiferromagnetic groundstate properties of the spin-1 fermionic system. This provides 
further insight in understanding the signature of the strongly interacting Bose- Fermi mixture 
and the two-component Fermi- Fermi mixture with population imbalance 32| . In this section, 
we use the TBA approach to investigate the groundstate properties and magnetism of the 
ID mixture. Application of an external magnetic field to the polarized fermions causes 
Zeeman splitting of the spin-up and spin-down fermions into different energy levels. The 
groundstate can only accommodate fermions that are in the lower energy level. Therefore 
it is expected that when the direction of the magnetic field is along the spin-up (H > 0) 
direction, spin-down fermions can no longer populate the groundstate. We shall verify this 
picture from the TBA formalism. Under a magnetic field, the system possesses three phases: 
I) a pure boson phase; II) a mixed boson-fermion phase; III) a fully polarized fermion phase. 

The derivation of the TBA equations for the model (TjQ) is standard [69|. Here we present a 
proper set of dressed energy equations which are convenient in the analysis of the groundstate 
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properties. These equations are 



e( * ) = *-*-2* + ^y_ fl ?T^ (33) 

0(A) = F, 1 ^ Ce " (fc) 



27r y-g ? + (A - fc) 2 
J_ /"* c^-(A) _ J_ r B 2c0-(AQ 
2tt y_ d ^ + (A-A) 2 2tt y_ B c 2 + (A - A) 2 1 J 

H ^ f B c0"(A) 

?/>(A) = -u b + — / -= — — . 35 

1 ^ ^ 2 27ri_ B f + (A -A) 2 1 ' 

The chemical potentials p, and /x^ are for the total number of particles and for bosons. The 



69| byexp(e(fc)/T) 



magnetic field is again denoted by H. The dressed energies are defined 
p h {k)/p{k), exp(0(A)/T) = a h {\)/a{\) and exp(^(A)/T) = r h (A)/r(A). Here p(k) {p h {k)), 
<t(A) (a h (X)) and r(A) (r ft (A)) are particle densities (hole densities) in k space, A space and 



69]. 



A space, determined by the Bethe ansatz equations (J6J) in the thermodynamic limit 
The superscript — denotes the negative part of the dressed energies, which correspond to 
the occupied states. The positive part of the dressed energies correspond to the unoccupied 
states. The pressure of the system can be obtained from 

1 f Q 

P = — / e~(k)dk. (36) 

27r J-Q 

The dressed energy equations (133]) - (135]) provide an elegant way to analyze the groundstate 
properties in terms of external fields. For H = 0, we see that e(k) has finite Fermi points. 
However, 0(A) and ip(A) have Fermi surfaces at infinity. After taking Fourier transforms in 
equations ( |34l) and ( |35l) we obtain the single dressed energy equation 

„n ,2 1 f Q 2ce~(k) 



42 



70| with 



which is equivalent to the dressed energy equation for the spinless Bose gas 
the chemical potential \xb — A 4 + £*&■ Therefore, the groundstate for the mixture with H = 
is exactly the same as that for Lieb-Lininger bosons. For strong coupling, the groundstate 
energy 

E Ion/ 4 12\ 

can be derived from (|37|) . We have thus shown that the groundstate at H = is populated 
only by bosons. 
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As the external field is turned on, the Fermi surface for 0(A) rises gradually. In this case, 
we have 



«» = *->-k H+ H\j$$r (39) 

« A) = _ w+ l Jr+ ££ ? ^_. (40) 

Analysis of these dressed energy equations reveals that pure Zeeman splitting between the 
two- component fermions does not energetically favour the fermions with spin-down. There- 
fore, the pure Zeeman field H can only drive the system from the pure boson phase into 
a mixture of fermions with spin-up and bosons or a phase of fully polarized fermions for a 
large enough external field. The mixture of partially polarized fermions and bosons requires 
two Zeeman splitting parameters which can maintain both fermions with spin-up and spin- 
down in certain regimes (see, e.g., the various regimes for ID three-component interacting 
fermions 711] ). As the magnetic field is increased, more fermions with spin-up populate the 
groundstate. Beyond a critical field value, the system will be entirely occupied by fermions of 
one species. A phase transition is expected to occur when the field exceeds the critical value 
Hff . When H > H^ 4 , 0(A) will become non-negative. From the dressed energy equations 
(159|) and (HQ)), we have 

w—k + I 11 (41) 

Now when H = H^ 1 , the Fermi sea for the spin pseudomomenta vanishes, i.e., 0(0) = and 
the density of bosons m& = M^/L = 0. Hence, the critical value is H^ 1 = 8P/c, where P is 
the pressure per unit length of the system, which will be determined below. 

For the mixture of fermions and bosons, the critical field is twice the value of that 
for spin- 1 fermions with a repulsive interaction discussed in the last section. The reason 
follows from the difference in their statistics. For the repulsive fermions, the groundstate is 
antiferromagnetic, with an equal number of spin-up and spin-down fermions. As the external 
magnetic field is increased it reaches the critical value H^, at which half the total number 
of fermions with spin-down are flipped. Whereas in the mixture of polarized fermions and 
bosons the groundstate is populated only by bosons. When the field exceeds the critical 
value > all of the bosons are driven out of the groundstate while all of the fermions 
are polarized. Therefore, full polarization in the mixture costs twice the energy of fully 
polarizing the spin-| fermions. 
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The dressed energy equations (139|) and (|40l) can be dealt with analytically for an external 
field < H < H c and in the strong coupling limit 7 ^> 1. Substitution of equation fj4*0|) 
into equation (|3T?j) and some lengthy iteration calculations with the relations e{Q) = and 
4>{B) = gives the leading behavior of the pressure per unit length as 



P 



2fi 3 / 2 HQ 1 / 2 
8Q 2 



tan 



2vr 

-1 



2 1 [2B\ 

1 -; tan (t) 



2B 

c 



+ 



3tt 3 c \ c J 3tt 3 (c 2 + AB 2 ) 

where we have adopted the notation 

2fi b 1 (2B\ 
n = n+ -^tan" 1 ( — ) • 

In the above equations the integration boundary for <p(\) is 



(42) 



(43) 



B _ j8P/c-H + 2fi b 



H-2fM b 



(44) 



and the Fermi point Q can be expessed as 

Q « fi ] 

2P 





1 - - tan" 1 










7T 





tan" 



1 /' 2 ^ 



+ 



APB 



fivrc \ c J ' nVL(c 2 + 45 2 ) 

The density of bosons m b can be estimated from the pressure 
m b = dP/dfib as 



« — tan 

n 7r 



¥) 



(45) 

2]) and the relation 
(46) 



Using the relations (jHJ) and (j46]) we have fi b « — APc~ l cos 2 { 11 ^-). The chemical 
potentials for bosons and fermions are chosen as /!# = ^ + fi b and /ii? = /i. Furthermore, 
on neglecting terms of order 0(l/7 2 ) and using the particle density relation n = dP/dfi, we 
obtain the pressure per unit length and the free energy in the form 

2 



P 
F 



-n 3 ir 2 



1 6 /m b | singg 
7 \ n ix 



, - , m b | sin(y 
7 \ n 7r 



2 3 

from which we can examine the magnetic properties of the model. 



(47) 
(48) 



14 



A. Magnetic properties 



We first see that the energy per unit length E = F + m z H is consistent with the results 
obtained in Ref. 51J. It follows that 

(49) 



H « — ttV 

37 



1 + cos (—) 



The susceptibility x follows as 

dm z 37 
dH 1671-% sin 



X(H) = = , a ^±, mb «v (50) 



n 



From the relation (HDJ, we find sin(^) = 2 (jjg-r) 1/2 (l - if&) 1/2 - As a result the 
susceptibility is given by 



where the critical field if^ 1 ~ ^Ep. It is clear to see that the susceptibility is divergent 

rM 



in the vicinities of if = and if = if* . The magnetization in the vicinity of these points 
belongs to the universality class of square-root field-dependent magnetization. This indicates 
a van Hove type singularity of the empty band. It is evident that pure Zeeman splitting for 
the mixture of polarized fermions and bosons may trigger three phases: a pure boson phase 
when the external field is absent; a fully polarized phase when the external field exceeds 
the critical value ff^ ; and a coexisting phase of fully polarized fermions and bosons for 
< H < H^' 1 , as shown in figure [31 

This phase diagram is reminiscent of that of ID weakly attractive fermions with popula- 



tion imbalance 



27 



72| . where the fully paired phase occurs only for the external field H = 0; 
the fermions are fully polarized for H > H c = n 2 (7r 2 +2|^|); the paired and unpaired fermions 
coexist for < H < H c . However, there are subtle differences with the phase diagram for 



ID strongly attractive fermions with population imbalance [19] . where the bound pairs in 
the homogeneous system form a singlet groundstate when the external field H < H c \- A 
completely ferromagnetic phase without pairing occurs when the external field H > H c2 
and the paired and unpaired atoms coexist for an intermediate field if c i < H < H C 2. The 
essential differences between the magnetic properties of the interacting Bose-Fermi mixture 
and interacting fermions with polarization are due to their different statistical signatures. 

Plots of the magnetization and susceptibility versus the external magnetic field H are 
shown in figures H] and [5l The susceptibility is infinitely divergent in the vicinities of if = 

15 



and H = H^ 1 . We can also see that as the interaction strength 7 increases the critical field 
tends to zero. If the interaction strength is weaker, the particles have more freedom to move 
along the line. Consequently, there will be a strong spin fluctuation which makes it harder 
to fully polarize the entire system. On the other hand, for strong repulsion the spins are 
"frozen" , thus making it easier to fully polarize the system. 

For the weak coupling regime, we obtain the groundstate energy of the mixture of fully- 



polarized fermions and bosons from the Bethe ansatz equations (jSJ), with result 



52] 



E 

T 



Mf M*c 2M b M lC 



(52) 



The magnetization m z « -^(y/2H + ^) and the susceptibility \ ~ sWh ^°^ ow fr° m 
this equation. These results show that in the weak coupling regime the square-root field- 
dependent behavior of magnetization emerges for finite external field. This is different from 
the antiferromagnetic behavior of the ID Fermi gas with population imbalance. 



V. CONCLUSION 

We have studied external field- dependent magnetic properties and statistical effects in 
ID polarized strongly repulsive fermions and in a ID mixture of polarized fermions and 
bosons with a repulsive interaction. We have found that the linear field-dependent behavior 
of the magnetization in the polarized Fermi gas persists for weak and finitely strong coupling 
regimes. The susceptibility is found to be finite. However, in the extreme limit 7 — > 00, the 
susceptibility is infinitely divergent due to its paramagnetic signature. The spin-spin inter- 
action is effectively described by the isotropic Heisenberg spin chain with antiferromagnetic 
coupling constant J ~ —AEp/^f. A quantum phase transition from the partially-polarized 
phase into the fully-polarized phase occurs when the external field is greater than the crit- 
ical value « 8Ep/ / -f, recall figure EJ For the weak coupling regime, the critical field is 
~ n 2 (ii 2 — 27). From these configurations, we can predict the subtle segments for the 
ID Fermi gas in a harmonic trapping potential: partially-polarized fermions lie in the center 



of the cloud whereas fully-polarized fermions sit in the two outer wings (see also Ref. [36| ) . 
The model of polarized fermions with repulsive interaction provides a tunable many-body 
system exhibiting novel critical behavior. It is highly desirable to probe this many-body 
physics through experiments with ID interacting fermions. 
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For the mixture of polarized fermions and bosons we have shown that the groundstate 
is only populated by bosons in the absence of a magnetic field. When the external field is 
applied, spin-up fermions populate the groundstate. This leads to an infinitely divergent 
susceptibility. The fully-polarized fermions and bosons coexist in the range < H < 
H^ 1 pa 16Ep/j. Another phase transition from this partially-polarized phase into the fully- 
polarized Fermi phase takes place when the field H > H^ 1 , recall figure [3j The susceptibility 
diverges. This phase diagram is somewhat reminiscent of weakly interacting attractive 



fermions 



271 ]. where weakly interacting BCS pairs can be viewed as strongly repulsive bosons 



in the mixture. This signature was recently observed in experiment 32j. There are subtle but 
essential differences between the mixture and the attractive fermions due to the different 
statistical signatures of the boson and the bosonic dimer. The mixture in a harmonic 
trapping potential has distinct segments: a boson-fermion mixed phase in the center of 
the cloud and a fully-polarized fermion phase in the two outer wings (see also Ref. |51|). 
These exotic magnetic properties may also possibly be observed in experiment with ultracold 
fermionic and bosonic atoms through photoemission spectroscopy techniques [73]. 

This work has been supported by the Australian Research Council. The authors thank 
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FIG. 1: Magnetization m z (H) vs external field H for n = 1 and different values of interaction 
strength 7. The small field magnetization curves are plotted from equations dH) and (|2ip for the 
weak and finitely strong coupling regimes, respectively. 
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FIG. 2: Phase diagram for ID polarized fermions with strong coupling c = 10. 
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FIG. 3: Phase diagram for the ID mixture of polarized fermions and bosons with strong coupling 
c = 10. 




magnetic field H 

FIG. 4: Magnetization m z (H) vs external field H for the Bose- Fermi mixture with n = 1 and 
different values of 7. The magnetization curves are plotted from equation (|49p . 
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FIG. 5: Susceptibility x(-ff) vs external field H for the Bose-Fermi mixture with n = 1 and different 
values of 7. The susceptibility is evaluated from equation ([50]) . 
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